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Qli| Abstract 

We study in detail a supersymmetric Peccei-Quinn model, which has a DFSZ and 
a KSVZ version. The fields breaking the Peccei-Quinn (PQ) symmetry correspond 
to flat directions (flaton fields) and have unsuppressed couplings when PQ symmetry 
^ ■ is unbroken. The models have interesting particle physics phenomenology. The PQ 

scale is naturally generated through radiative corrections; also, in the DFSZ case, the 
\x problem can be solved and neutrino masses can be generated. Cosmologically they 
lead to a short period of thermal inflation making the axion an excellent dark matter 
candidate if one of the flaton fields has a positive effective mass-squared at early times 
but with too low a density in the opposite case. A highly relativistic population of 
axions is produced by flaton decay during the subsequent reheating, whose density 
is constrained by nucleosynthesis. We compute all of the relevant reaction rates and 
evaluate the nucleosynthesis constraint. We find that the KSVZ model is practically 
ruled out, while the DFSZ model has a sizable allowed region of parameter space. 



1 Introduction 



With the discovery of the instantons it was realized that the pure gradient topological term 
of the QCD Lagrangian 9qcd9 2 s/32'k 2 FF can generate important physical consequences. In 
fact the induced CP violation affects the electric dipole moment of the neutron suggesting the 
limit 9 qcd < 1CT 10 . The most attractive explanation for the origin of such a small parameter 
would be the Peccei-Quinn mechanism |TJ. There is supposed to be a spontaneously broken 
global £7(1) symmetry (PQ symmetry), which is also explicitly broken by the color anomaly. 
The corresponding pseudo-Goldstone boson is called the axion 0. 
The PQ symmetry acts on some set of fields <pi with charges Qi, 

0._> e «fc«^. (1) 

To generate the require U{1)pq x SU(3) c x SU(3) c anomaly we must choose appropriate 
particle spectra. Depending on the charge of the SM matter content, we can have the KSVZ 
(hadronic) models || in which only some extra heavy quark fields are PQ charged or the 
DFSZ H models in which, beyond the extra matter content (at least two Higgs fields), also 
the Standard Model (SM) matter is PQ charged. 

Denoting the vacuum expectation values of the scalar fields charged under PQ by t>i/\/2, 
the PQ symmetry breaking scale F PQ is defined as Fp Q = J2iQl v f, an d the axion mass is 
given by FpQin a = (79MeV) 2 iV where N is the number of quarks with PQ charge. Col- 
lider and astrophysics constraints require F P q > 10 9 GeV, and defining as usual M P1 = 
(SitGn)- 1 / 2 = 2.4 x 10 18 GeV this allows the range 

10 9 GeV<F PQ <M P1 . (2) 

With typical assumptions about the cosmology, the requirement that axions give at most 
critical density places Fpq towards the bottom of this range. In the particular case that the 
axions are radiated by strings with no subsequent entropy production, one probably requires 
@ F PQ ~10 10 GeV. 

In a model with unbroken supersymmetry, the holomorphy of the superpotential ensures 
that PQ symmetry Eq. (|I|) is accompanied by a symmetry acting on the radial parts of the 
PQ charged fields 

fa^e^fa. (3) 

The corresponding pseudo-Goldstone boson is called the saxion (or saxino), and the spin-half 
partner is called the axino. Soft supersymmetry breaking gives the saxion a mass of order 
100 GeV, and the axino typically has a mass of the same order || though it may be very 
light in special cases.Q 

Let us consider the potential of the fields fa which break PQ symmetry. In a supersym- 
metric model there have to be at least two, but let us pretend for the moment that there is 

1 We assume gravity-mediated supersymmetry breaking, which typically gives soft scalar masses in the 
range 100 to 1000 GeV, and use the former estimate for definiteness. 
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only one. Its potential will be of the form 



1 00 Ul 2n+4 

V = V - m 2 |0| 2 + V + £ A n ^L_ . (4 ) 

^ n =l JW P1 

The non-renormalizable terms are expected to have coefficients A n ~ 1. If the renormalizable 
coupling A is also of order 1, m ~ -Fpq and the non-renormalizable terms are negligible. 
(Remembering that there are at least two complex fields, only the particular combination of 
fields corresponding to the saxion will have the soft mass of order 100 GeV.) In this sort of 
model one can hope to understand a value Fpq ~ 10 10 GeV since that is the supersymmetry 
breaking scale [0, but it may be hard to understand a bigger scale. 

We are concerned with a different class of models [[8], [||, [10], 11, 12], in which </> represents 



a flat direction of supersymmetry. The quartic term is then absent, while m is a soft 
mass of order 100 GeV. If the 6 term is present with unsuppressed coefficient the vev 
is (0) ~ A x ^ 4 10 10 GeV. If instead the 8 dominates one has ((f)) ~ A 2 1 ^ 6 10 13 GeV and so 
on. For future reference, note that the mass of \<f>\ in the vacuum is also of order of m, and 
that the height of the potential is given by 

T/l/4 \ / / ,\ \ 1/2 

\ - f ^ \ . (5) 



10 6 GeVl V 10l ° GeV / ' 

Fields of this kind, characterized by a large vev and a flat potential are called flaton fields, 
and the particles corresponding to them are called flatons [[0| . 



In a supersymmetric model there are n > 2 complex fields which all acquire vevs. We 
are interested in the case where these are flaton fields. Then there are 2n — 1 flaton particles 
with mass of order 100 GeV, and n flatinos with typically similar masses. The saxion (axino) 
is a linear combination of the flaton (flatino) fields, with no particular significance. 

The rest of this paper is as follows. In the next section we define our models and 
summarize the cosmology. In section ^|, we give the general structure of the flaton and 
flatino masses. In section || we analyze the general self interactions between flatons and 
flatinos. In section |^ we see the effect of the interaction of the flatons with the matter fields 
(The KSVZ flatons interact only with gluons and gluinos whereas DFSZ flatons interact also 
with ordinary matter and supermatter.) In section ^| we find the parameter space regions 
that can satisfy the cosmological constraints. We conclude in Section |7]. 



2 The model and its cosmology 

We consider a DFSZ and a KSVZ (hadronic) model. 

In both models, there are two flaton fields P and Q, interacting with the superpotential 



11, 12 



W &aton = ^-P 3 Q. (6) 
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In the hadronic version, the interaction with matter is 



Wflaton-matter — ^EiEiE\P (7) 

where Ei and Ef are additional heavy quark and antiquark superfields. 
In the DFSZ version, the interaction is 

WW_ matter = IXNNP + yt-E\EiPQ (8) 

where N are the right handed neutrino superfields and H 12 the two Higgs doublets. Due to 



the second term we can provide a solution to the \i problem |L4| . In such case we can add to 
the superpotential of the minimal supersymmetric standard model also the terms h u lH 2 N 
that generate the necessary mixing between left and right neutrinos to implement a see saw 
mechanism which can explain the solar and atmospheric neutrino deficits. 



The cosmology of the DFSZ model has already been considered in a rough way ||15|| . Here 
we relatively complete treatment of both models. 

To study the cosmology of the flaton fields, we can safety analyze only the superpotential 
Wflaton for both models. With the inclusion of the soft susy breaking terms, the potential is 



V = ml\M 2 + ™|I0q| 2 + ( 9 I<M 4 |<^| 2 + I0p| 6 ) + (j^f^Q + h - c ) • ( 9 ) 

The soft parameters mp, rriQ and Af are all of order 10 2-3 GeV in magnitude. It is assumed 
that mp and w?q are both positive at the Planck scale. The unsuppressed interactions of <pp 
give radiative corrections which drive mp to a negative value at the PQ scale, triggering a 
vev for (pp. As a result, when 0p gets a vev, the tadpole term proportional to Af generates 
automatically a vev for the 4>q field, both of them are v q ~ vp ~ 10 10-12 GeV. 

In the early Universe when H > 100 GeV, there will be effective values mp(t) and m^{t). 
Supergravity interactions will make both of them at least of order if 2 in magnitude.^ The 
cosmology depends on the signs of these effective values. If mp(t) is positive, <p P is held at the 
origin by the finite-temperature potential until T ~ |j7ip|. But the potential Vq dominates 
the energy density in the regime |m P | <T< V^ 4 , leading to about ~ In (K) 1//4 /|mp|) ~ 10 



e-folds of thermal inflation |13| 



To discuss what happens after thermal inflation, we suppose first that Triqit) is also 
positive in the early Universe, so that 0q is also trapped at the origin during thermal 
inflation. When thermal inflation ends, <pp moves away from the origin, which destabilizes 
cf)Q. The fields (ftp and 0q move around an orbit in field space, which would be closed if there 
were no energy loss. If the only energy loss came from Hubble damping the fields would 
oscillate back and forth many times around an almost-closed orbit. However, the parameter 
determining the strength of parametric resonance is q ~ g^o/m ~ 10 8 , where g ~ 1 is a 
typical coupling, <3> ~ 10 10 GeV is the amplitude of the oscillation and m ~ 10 2 GeV is its 
angular frequency. One therefore expects that parametric resonance will efficiently damp 



"During inflation this result might be avoided (say by D-term inflation) but it should still hold afterwards. 
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the orbit, converting most of the energy into flaton particles. At first these particles are 
marginally relativistic, but after a few Hubble times they become non-relativistic. The rest 
of the energy resides in the homogeneous oscillating flaton fields, now with small amplitude 
and therefore almost simple harmonic motion corresponding to some more non-relativistic 
flatons. The flatons decay leading to final reheating at a temperature | 13| 

10 n GeV\ / m ^ 3 



where m is the mass of the lightest flaton and g^n ~ 100 is the effective number of relativistic 
species at reheat. 

Two different axion populations are produced. One population is radiated by the PQ 
strings that form after thermal inflation. They become dark matter with abundance H 

~ (-^pq/10 10 GeV) 1 ' 2 . The present scenario predicts Fpq ~ 10 10 GeV, making the axion 
an excellent dark matter candidate. Note that in contrast with the general case, the axion 
density in this scenario cannot be reduced by entropy production after the epoch T ~ 1 GeV 
when the axions acquire mass; pre-existing long-lived particles that might do the job have 
been diluted away by the thermal inflation. 

The other axion population, that is our main concern, comes from the decay of the flatons 



13], M. This population is still relativistic at nucleosynthesis, and its density must satisfy 



the constraint 

< 5N U ~ 0.1 - 1.5, (10) 




NT 

where p u denotes the energy density of a single species of relativistic neutrino and SN U the 
number of extra neutrino species allowed by nucleosynthesis. 

If there were just one species of flaton field <f), this would give the bound [[HJ 

1/3 




NT 



f 


^43/4 




\9rh 


f 


f43/4 




\9rh 



B a 



l-B a 

1/3 

B a , (11) 

where B a = T a (0 — > a + a) /T tat is the branching ratio, and in the last line we assumed 
B„ C 1. From Eq. (|ToD we get the bound 

T(4,^X) 1.5) {43/4 ) ■ {U] 

Varying Tr# from 6 MeV torn ~ 100 GeV we get a factor two variation coming from the 
number of degrees of freedom. For 8N U varying from 0.1 to 1.5 we get B a < 1/3 to B a < 0.02. 
Our model has three flaton particles and the quantity B a to be used in Eq. ( |Tl"D is 

m 



B a = y~] rjBj where rj 



r> r(I^a + a) + 1 1 T(I-,a + X) 

Bl = f^iT) • (13) 
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Here rij is the number density of the Jth flaton just before reheating, related to the mass 
density pi by rij = pi/mi. In principle one could calculate the nj (coming from parametric 
resonance and some residual homogeneous flaton oscillations decay) but we have not done 
that, and to estimate the allowed parameter range for the model we shall simply assume 
77 = 1/3. 

Next consider what happens after thermal inflation in the case that Tnq(t) is negative. 
During thermal inflation, while <pp is trapped at the origin, \4>q\ will have some value 3> -Fpq 
determined by a higher-order non-renormalizable term. At the end of thermal inflation, all 
three flaton particles will be produced in the manner we have described, and the nucleosyn- 
thesis constraint still holds. The difference from the previous case is that axionic strings are 
not produced, so that dark matter axions are produced only by the quantum fluctuation of 
the axion field during inflation. The density is now |J, §fl a ~ 3(#/7r) 2 (F PQ /10 12 GeV) 12 , 
where 9 < it is the misalignment angle. Again, thermal inflation means that the abundance 
will not be diluted by entropy production, so discounting an accidentally small 9 we again 
have a rather definite prediction, which is too low. For the axions to be the dark mat- 
ter in this scenario, one would have to increase Fpq by generating it from a higher-order 
non-renormalizable term. 

Finally, if m 2 P (t) is negative in the early Universe, there is no thermal inflation and one 
is back with all the uncertainties of more general models [|Tj], [18|, [L9[ . Axion cosmology now 



depends on the scale of the inflaton potential, the reheat temperature after inflation, the 
decoupling temperature of the flatons and so on. We have nothing to say about that case. 



3 Flaton and flatino spectrum 

3.1 Flaton spectrum 

We write the flaton fields as 

Vp + P i^P 

-e v p 



V2 
vq + Q 
V2 



Ar 



e "o , (14) 



and we shall take vp, Vq, Af and / as the independent parameters in the potential Eq. (pj). 
The main components of the axion field are 

VP A P + 3^A Q (15) 



Fpq F P q 

where Fpn = vp + 9vq and we have neglected the O (vew/Fpq) components along the 
directions. The orthogonal field to the axion (both are CP Odd) corresponds to a flaton 
particle. It is 

* = -£-AQ-3g-Ar (16) 

-TPQ -TPQ 
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which has a mass 



where 



Ml 



^ Af n V / FpQ = --fiAf (x 2 + 9 ) . (17) 
2M Pl v Q <T f y ) K ' 



V _ V P V Q ns x 

~g = (18) 

Up . , 

X = — . 19 

For future convenience we have introduced a quantity fi, related to the g appearing only in 
the DFSZ model. At this stage results depend only on the ratio fi/g defined by (0) and 
they apply to both models. Since M\, is positive, Af and / must have opposite signs. 

The other two flaton particles correspond to the CP even fields P and Q. They have a 
2®2 mass matrix whose components are 



2 



M QQ = (2°) 

Mpn = 9/^-3-^ = 
PQ J M&x 9 + x 2 

M|p = 3/ 2 p ;, 9 ; - 3 ^ = \T— (x 2 + 3) /x 2 - 3 ^ 

pp J M| x x 2 9 + x 2 # 2 V / P 9 + x 2 




Here two mass parameters m 2 P , rriQ in Eq. ([J) are replaced in favor of Vp, vq. Performing the 
rotation from the flavor basis \\P Q\\ to the mass basis \ \Fi FzW 

P = cos aF2 — sin aF\ 

Q = sin aF 2 + cos aF 1 , (21) 

the mixing angle results 

2M po x , , 

tan 2a = — = = -6— = -a . 22 

Note that tan 2a ranges from — v^3 for x = v^3 to zero for x ^> 1 and x — > 0. For each x 
there are two solutions given by 

cos 2 «i = L_) (23) 



2 I VT+ (T 



1 / . 1 



cos ' 2a2 = 2( 1 + 7TT7j (24) 

with sin«i 5 2 > 0. The two solutions will be relevant when we discuss the parameter space 
for the decay of flatons into axions. 
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The two eigenstates have masses 



M 2 2l = ^-( f -(12(x 2 + 3)^ + (3 - x 2 )^) ± \ f -{\2 f - + ^)|v / 9 + 42x 2 + x 4 ) (25) 
' 2 ^ g /I g g /I J 

The requirement of a positive definite spectrum {m 2 Fl > 0) gives the constraint 

qA f 9 + x 2 . . 

yi<y = ~TjT^ r<y ' 2 ' (26) 

where 

y 1;2 = (21 + ± V9 + 42x 2 + x 4 ) (27) 

or 

i (21 + x 2 - V9 + 42x 2 + x 4 ) < -^t- 1 < \ (21 + x 2 + V9 + 42x 2 + x 4 ) (28) 
and requiring positive diagonal elements implies also 



(x 2 + 3) (x 2 + 9) 
y < y 3 = ^ ^ ^ 29 

x^ 



3.2 Flatino spectrum 



From the superpotential Wflaton we can directly extract also the flatino's mass matrix whose 
eigenvalues are 

q iwf 2 f 2 

M| =^^[x 2 + 2±v / i T TT]=9^/i 2 [x 2 + 2± v / x^TT] (30) 
F2 ^ Ay x 2 9 

The eigenstates Fi, F 2 are related to the flavor states P, Q by 

iq = cos aP + sin ciQ 

F 2 = — sin aP + cos <S<5 (31) 

where tan 2a = x, 

A parameter space analysis indicates that we have always M Fl > 2M Fi . This automat- 
ically forbids the decay of F\ to flatinos leaving open only the decay into flatinos of the 
heavier F 2 and ip' flatons. 



4 Interactions between flatons 

Now that we know the general mass matrix structure of the flatons and flatinos, common 
both to the DFSZ and KSVZ models, we can start analyzing the various decay rate between 
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flatonic fields.^ We begin with the decay channels induced by the kinetic term and the 
superpotential Wflaton, which are common to the KSVZ and the DFSZ models. 



4.1 Derivative and cubic interaction terms between flatons 

The flaton interaction terms with at least one derivative are given by the lagrangian 

2v P P + P 2 . v% , 2 , n v\ 



(da) 2 + 9-5- {di)'f + G^^d^'da] + 



2^P -PpQ -PpQ -PpQ 

2 ,,2 „,2 



2V Q Q + Q 2 V-q 2 2 „VQVp / . 

9 „ 2 [9-dr(0a) + pa-W) -6-^2— 9^ 9a] (33) 

zy Q ^PQ ^PQ ^PQ 

From this expression we can extract the following terms expressed in mass eigenstates. 
The trilinear derivative interactions with no axions 

(dip') 2 [(-9 sin a + x 3 cos a) F x + (9 cos a + x 3 sin a) F 2 ] , (34) 

FpqxVx 2 + 9 v 7 v ' 

The trilinear derivative interactions with only one axion 

6 

Lf-^'o, — dtp' da [(coso; — xsina) F 2 — (sin a + xcosa) Fi] , (35) 

F PQ Vx 2 + 9 

The trilinear derivative interactions with two axions 

LF aa = \dfj,a\ 2 . ((9 cos a — x sin a) F% + (9 sin a + x cos a) F 2 ) . (36) 

FpQVr + 9 

All the above derivative interactions can be transformed in scalar interactions if we are 
working at tree level and with on-shell external particles 

0i (^0 2 ) = \ « - Ml - Ml) 0!0 2 03 (37) 

The cubic interactions come also from the flatonic superpotential and the soft terms 

T ( 9 f2 V P V Q | 5 2 VP M 2 ,VqX \ 3 | (27 f 2V>Q o M j' X \ p2 n 

4,3 " ^2 J M 2 ,2 J M 2 l v 2 P (x 2 + 9)J + \2 J M 2 X v P (x 2 + 9) J W 

+ k pA) pq 2 + (v^) pw + (LMkp) qw , (38) 

\2 J M 2 J V ^4M p1 ^qJ ^ \4 M Pl v 2 Q y ' 



3 The two body decay rate is given by the expression 

r(»-»/i/ 2 ) = J ^cos6A{Ei = £ 2 ,cos0) (32) 

where S=l if /i ^ / 2 and 5 = 1/2 is = f 2 ■ E x = E 2 = (a/. ( 2 - M* - M%) /2M % , \P\ = 
1/ (2Afi) y^Mf - (Af/, + M/ 2 ) 2 ) (m? - (A/ /x - M /2 ) 2 ) and A (..) is the square amplitude of the decay. 
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In the mass basis, we get 



Mm sin a 



(— 2x cos 2 a + x sin 2 a + cos a sin a ) + 



6 f 2 ii 2 \/x 2 + 9 / n o o 2 

5 — I —18 cos a sin a a; + 9 sin aa; + 3 cos a x 

g 2 xF PQ v 

- cos a sin 2 a(-9 + a; 2 )) } F 2 F 2 = — 



For the full trilinear F 2 i\)'i\)' interaction, we have to add up the terms in Eqs. 
to obtain 



(39) 
and PP 



Ml 



'F 2 i>'V - 7) /n i 2 1? 

2xV9 + x 2 F PQ 
(9 cos a + sin ax 3 ) ^1 - F ; 



(3 sin a + x cos a) (x 2 + 9) 



M%_ 



(40) 



Af *w^f 2 



Collecting the above formulae one finds the decay rates among flatons and axions 



1 M| 2 

n — 7 — n r [x cos a + 9 sm a ) 

32ttF 2 q (x 2 + 9) 1 ' 

1 Ml , 2 

1 '— xsina + 9 cos a) 



32vrF 2 Q (x 2 + 9) 
1 



Ml 



1-4— ^U, 12 



T (F 2 — > aa) = 
T (Fx — > aa) = 

T (F 2 -> FiFx) 

r (f 2 -> vV) 

r(F2 ^^ K 16^ 2 Q (x 2 + 9)r M^ 2 
T -> aF 2 ) = T (F 2 -> e#' 

T -> aFi) = T (Fx -> a^' 



1-4- 



32ttM F2 ^ M 2 2 
1 



- 2 —\^F 2 ip'i}'\ 



(41) 



1 — t-^t I (3 cos a — 3xsmaJ 



•0' *1 



Energy conservation will of course forbid some of these reactions, depending on the flaton 
masses. As M Fl < Mm the channels Fi — > if;' if)' and Fx — > ip'a are always forbidden. 



4.2 Interaction terms between flatons and flatinos 

The trilinear Lagrangian terms responsible for the decay of flatons or flatinos are 

= J-LpQ + , q P)PF + ^^^ 

^Mpi \ fipQ ^ P Q J 
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SJ '2Pv P PQ + 2i^- (3v Q if;' + vpa) F 7s Q ) (42) 



2M P1 V ^PQ 

(the tilded fields are the fermionic superpartner of the respective P and Q scalars). Let us 
denote the Yukawa couplings between the flaton (or the axion) and the flatinos in mass basis 
by —Lyuk = Yijk4>iFj (1, 75) -Ffc/2 where 75 is taken for fa = a, ip'. We find from Eq. fl42|) the 
following expressions for the Yukawa couplings 



6ffiy/x 2 + 9 . 2 _ . 

r p p 5 = ( x cos a — sin a ) cos a — x sin a sin 2a 

6ffiVx 2 + 9 . . 2 _ . 

r p 5 p = (x cos a — sin a) sin a + xsinasinza 

F lF2 F 2 gxFpQ LV J j 



6//xVx 2 + 9 r . ~ . 1 / . \ • q ~i 

j p p p = —x sin a cos 2a H — (sin a — xcosa) sin 2a 

FlFlF2 gxFp Q 1 2 v 1 J 

6/^^ + 9 . 2 _ 

r p p p = isina + cosa) cos a + x cos a sin 2a 

F2FlFl gxF PQ LV ; J 



6ffxVx 2 + 9 . . 2 _ . Q ~! 

V p 2 F 2 F 2 = ^ K X Sm a + COS a ) Sm a — X COS a Sm 



QffiVx 2 + 9 1 
*f p p = +xcosacos 2a (cosa + xsma) sm 2a 

F 2 F 1 F 2 gxFp Q 1 2 V ' J 

= —^—[— (3 + x 2 ) cos 2 a - 3x sin 2a] (43) 

6//1 / 2\ 2 

V ^2^2 = — [- (3 + x J sin a + 3x sin 2aJ 
y Vf x f 2 = gjp [-3xcos2a + ^ (3 + x 2 ) sin 2a] 
Yip p = — Ct— f2x cos 2 a — x 2 sin 2a] 
Y a p p = — [2x sin 2 a + x 2 sin 2a] 



6//i r > 

aFjifij — gxFpq 



Kp p = — - — [— x cos 2a — x sin 2a] 
From this we can extract the decay rates for Fj — >■ FjF k , if)' — > FjF fc , or F 2 — >■ FiF* (■*/>', a) 



r f f —> f f ^ — s (1 W*±¥*l Y(i {Mf ^ m ^ 2 Yy 2 

F (F t ^ F&) - —S \1 ^ j \1 j Y FzF]Fk 

m = *±8 (1 - m^-) 1 (1 - M 1 K 2 



■'■ r v= — m| — ; v — M 2 * — ) 
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1 1 
2\ 2 / / n/r _ i\ / r_\2\ 2 



16tt I M| /I M% ] \\ M* ) Ml ) F ^ F2 





\ Ml 








\ 2 M l> 







where 5 is a symmetric factor (1/2 for identical final states or otherwise 1). 



5 Interaction of flatons and flatinos with matter fields 

Now we study the interactions of the flatons with matter and supermatter, specified by 
Eq. fl7|) for the KSVZ case and by Eq. ([8]) for the DSFZ case. Through these interactions the 
over-produced flatons or flatinos could decay into ordinary matter while the number of the 
decay produced axions are sufficiently suppressed satisfying the nucleosynthesis limit (|10D. 



5.1 KSVZ model: Interactions between flatons and gluons 

We begin with the hadronic models in which the only decay mode available for the flatons is 
into two gluons coming from the anomaly (when the space phase will be available, we have 
to take into account also the decay into massive gluinos, in this discussion we neglect such 
a possibility). The respective one loop corrected decay rates are 



PQ 

and 



7T 



(45) 



T ( F ^9 + 9) = ^^N E -^-{ X +9)cosa[l + -——) (46) 

where Ne is the total number of the superheavy exotic quark fields (Me = h F vp 3> MfJ. 

We do not consider the flatino decay into a gluon and a gluino which will be irrelevant 
for our discussion. 



5.2 DFSZ model: Interactions between flatons/flatinos and ordi- 
nary matter 

The decay properties of the flatons in the DFSZ models involves the direct interactions 
between flatons and ordinary matter and supermatter In general the interaction between 
flatons and Higgs fields are quite interesting due to the fact that this two sectors, after the 
spontaneous breaking of the PQ and the EW symmetry, mix together. We notice that the 
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Peccei-Quinn symmetry prevents the introduction of a SUSY invariant mass term jj,HiH 2 , 
solving automatically the so called /i mass problem as mentioned before. 
Let us start by writing the Higgs-flaton potential 



V(H,, 



\Hi\ 



m 2 Hl + 



Mpi 

P0Q 



+ \H 2 



m k + 



+ 3r 0P 2 |0Q| 2 



Mpi 

2' 



+-A9 2 + g' 2 



pi 

2 



+ f 



y p 2 |0p 



+ c.c. 



(47) 



\H 2 



When the fields 0p,q get vevs, the m 2 HiH 2 mass term is generated dynamically. The size 
of such a term is fixed by 



wij = /i \ A g + —fi (x 2 + 3 



(4£ 



In the limit |m§| ^> the masses of the pseudoscalar A , of the CP even scalar Higgs field 
H° and of the charged Higgs fields H ± are almost degenerate 



m A°,H°,H ± 



sin (3 cos /3 

so from the constraint of positivity of such a masses we get 



^ + ^(x 2 + 3) < 



9 



(49) 



(50) 



In such a limit we also know that the mass eigenstate of the CP even electroweak sector 
H°, h° and of the CP odd one A , G° are 



H° = - sin (3 hi + cos (3 h\ 

h° = cos (3 hi + sin/3 h° 2 

A = sin/5 A° + cos/3 A° 2 

G° = cos/3A°-sin/3A° 



(51) 



where H\ — 4= (t^ + + ^4?) and H 2 = {y 2 + h 2 + iA 2 ) are the gauge eigenstates and 
tan/? = v 2 jv\. To allow the flaton decay into A , we want it to be light so that small tan/3 
is preferred in our discussion. Hereafter we will take tan/3 = 1. 
From Eq. (fl8D , we find 



Agfl 



+c.c. 



2' 
P 



!/.■/,/■ = ~f [K + K +Af + Af)[i- + ^ 

Q .x 2 + 3 lf 



Vp Vq 



xF\ 



PQ 



/ 9 \ v p 



P Q_ 

Cl> Vq 

V Q 



h°X - A" A") + 1 [ti{A% + h%Al)^2) 



|0 /i0 



,0 A o 



,0 A 



xFi 



PQ 



+ i—^—if/ ) + x 2 -fx 2 
9 



'P .6 
4 h 1 

v v P xFpq 
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It is then simple manner to get the decay rates for the kinematically more favorable decay 
channels F 1)2 -»• h°h° and t/j -> /*°^° 

r - ^ _ ^ (* 2 + 9) I* (, 4i^oV /2 , |2 

r ^ 2 ^ ^ j _ ________ ^1 _ _j (53) 

= ^5_4 f 1 - ___^__r fi _ _^_gy /2 w 



where 



— + 6— ) (x cos a — sin a) — Ax 2 — 



Ap^hh = sin2/3[ — + 6— (xcosa — sin a) — Ax — sin a] + 2 (xcosa — sin a) 
\ V 9) 9 

A-Fohh = sin2/3[f — + 6— ) (a; sin a + cos a) + Ax 2 — cos a] + 2 (2; sin a + cos a) 
\ f* 9 J 9 

^„=f_-a^M_ (54) 

If flatons produce a large number of flatinos, flatino decay into axions has to has to be 
suppressed as well. Primary importance is the production of the lightest flatino F\ which 
cannot decay into other flatons (axions) or flatinos. The flatino decay into ordinary particles 
comes from the superpotential W = -^H\HiPQ. We find that the flatino decay into a 
Higgs and a Higgsino (more precisely, the lightest neutralino Xi) nas the rate; 

1/2 



h° 



where = (sina + coso;)A^ xl and C^ 2 = (cos a — x" 1 sin a)N Xl . Here JV X1 denotes the 
fraction of lightest neutralino in Higgsinos. 

Let us now consider the flaton decay into ordinary fermions or sfermions. The mixing 
terms between flaton and Higgs fields allow a direct tree level coupling (after full mass matrix 
diagonalization) between the usual fermions and flatons. Parameterizing such a mixing with 
the parameter Q FH the effective flaton-fermion interaction is hfOpn so that the rate of decay 
is 

where N c is a color factor for the fermion /. Since $fh — (l^)> f° r light fermions 

(2m f < M F ) r (f* -> / + /) /r (Fi^a + a) ~ h 2 f v 2 EW /M 2 F . ~ m 2 f /M\ which is less than 
one. Therefore, the rate of the flaton decay into ordinary fermion cannot be made sufficiently 
larger than that into axions. 
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For the coupling between sfermions and flatons, we have two contributions. One is a 
direct coupling coming from the scalar potential 

V Fff = J^ ^ 2 + 9 ^L (h d tfm/3D L *D R * + h e t^(3E L *E R + h u U L *U R ) (56) 

(F% (x cos a — sin a) + F 2 (cos a + x sin a) ) + h.c. 

where D* denote down-type squark, etc. 

The other arises from an indirect coupling induced by the mixing between Higgs and 
flaton fields as for the fermion case. Taking in consideration the cubic soft A-terms we find 

V eff = h d A d 9 Fi Hl FiD L D R + h e A e 9 Fi Hl FiE L E R + h u A u 9 Fi h 2 FiU L U R (57) 

so that effectively we have couplings of the size 

Gfs} ermion ~ hf(fl + Af)-^ (58) 

-^PQ 



Diagonalizing the sfermion mass matrix we can write f R fi = On/i/i + ^22/2/2 + a i2/i/2 
(where an oc hf so for hf — ► we have a 12 — > 1). Considering the decay of the light flaton 
we get 



Q2 

T (F 1 ^fi+ 1) ~ N c - F l a 2 



2 

Tfl ~ 

■ . 1 - (59) 

As observed in Ref. [ I5| , the flaton may decay efficiently to two light stops as h t ~ 1 and 
~ 1 and thus a large splitting between light and heavy stops helps increasing the flaton 
decay rate to light stops. This kind of mass splitting occurs also in the Higgs sector and 
furthermore the light Higgs (h°) is usually substantially lighter than the heavy Higgs (H°) in 
the minimal supersymmetric standard model. This should be contrasted to the case with the 
mass splitting for stops which requires some adjustment in soft parameters. In this paper, 
therefore, we will concentrate on the flaton decay into Higg SGS clS db dominant mode of flaton 
decays. 



6 Parameter space analysis 

Our task is now to find the parameter space for which B a gets small enough. As a reference, 
let us try to see if B a < 0.1 (corresponding to 5N U < 0.7) can be obtained imposing the 
stronger condition that Bj < 0.1 for each / = Fi^,ip'- We first note that the decay rates 
calculated in the previous sections are functions of the 4 variables x = vp/vq, f /g, Af/fj, 
and Ag/fj, disregarding their overall dependence on F P q. 

Before starting our discussion, let us make some remarks. We are dealing with two kind 
of PQ models with a natural intermediate scale 1) the DFSZ and 2) the Hadronic model 
(KSVZ). 
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i) They have a common flaton potential and thus the same flaton and flatino spectrum. 
But they have different interactions between ffatons and matter. 

ii) The symmetries and parameter space constraints impose that the following decays 
are forbidden ip' — > aa, Fx — > ip'a, Fi — > ip'ip', -^i ~ > FF. 

Neglecting for the moment the model dependent flaton-matter interactions, we have to 
analyze the decays F 1]2 — > a a, F 2 — > ip'a and the orthogonal ip' — *> F\$ a plus the flatons- 
flatinos interactions that through the chain Flaton — > Flatinos — > Flatino — axion can 
also generate a non negligible axionic density at nucleosynthesis time. 

The decays of the flatons into two axions doesn't have any phase space suppression. If we 
choose a — a% (see below Eq. ^) then we can suppress the rate T {F\ — > a + a) only taking 

jrf 3 

x 2 = 18 and cosai|a.2 =18 = 0.426 (giving for example T (F 2 — > a + a) = 3.6 32n p 2 )■ On the 

PQ 

other hand, if we choose a = a 2 we can suppress only T (F 2 — > a + a) in the region x 2 = 18 

M 3 

and cosa 2 L.2 =18 = —0.905 (giving for example T (F 2 —> a + a) = 3.6 „ 9 p 2 )■ 

oz7r_TpQ 

The decay rate into a single axion has the phase space suppression constraint m 2 p 2 > Ml, 

that translated in our parameters reads — j^j- < 12. 

When the flaton's branching ratio to flatinos becomes sizable we have also to impose 

B (i -> F 2 ) B (P 2 -> aF 1 ) < 0.1 . (60) 



with B(l^F 2 )=2B (F 2 -> F 2 F 2 )+B (f 2 -> F 2 Fi) +2B(#^ F 2 F 2 ) +B -> F 2 Fi 

and also the important requirement that F\ has to be heavier than the lightest neutralino 
(Xi) and the light Higgs (h°) since it has the unique decay mode to a neutralino and a Higgs. 
Assuming mp l is much larger than the light Higgs mass, we impose a strong condition 

mp > fi that translates into the approximate bound 



6.1 Parameter space of KSVZ models 

As discussed already, the light flaton Fi can decay only into two gluons, so the simple 
requirement Bp — T[Fi — > gg\/Y[F\ — > aa] ^> 1 impose a quite strong constraint on the 
parameter space. This ratio is given by 

TjF^g + g) = qj (M Fl ) 4^ 2 / S ma(x 2 + 9) \ 2 / + 95 a s {M Fl ) \ _ 
r (Fi — > a + a) tx 2 9 E \x (9 cos a — x sin a) J \ 4 n J 

In order to have a large number we can use only the x parameter and the number of exotic 
quark ( to be as large as possible). For as (Mp) — 0.1, this ratio can be larger than 10 
accepting a tuning of the x parameter as follows 

- 0.04 < x < VlS + 0.04 for N E = 1, 
Vl8 - 0.3 < x < Vl8 + 0.3 for N E = 9. (62) 

However, in such region we don't find any solution in which, B a (F 2 ) and B a (ip') are less 
than 0.1 at same time. Therefore the KSVZ model cannot give a satisfactory solution of the 
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nucleosynthesis bound unless some extra fine tuning on the initial densities of the flatons 
F 2 , if)' can be dynamically achieved. 



6.2 Parameter space of DFSZ models 

In the DFSZ model the number of the possible decay channels is much larger than the 
other model and it is not an easy task to find, in the four parameter space, some easy 
understandable available region. To be as independent as possible of the soft supersymmetry 
breaking parameters we will try to make analytic computations on the rates of the flaton 
decays into Higgs particles, in particular into the lightest Higgs (h°) whose mass has an 
automatic upper bound of ~ 140 GeV p0[. We will concentrate on the region with - 



negative and |^ <C 1 and x ^> 1 which turns out to be required for B a < 0.1. 

To open the decay channels of the flatons into Higgs particles we have, in particular, to 
impose > Ma > which requires (from now on we will use the unequal symbols as 



strictly satisfied) 



< 



x 



with 
9 V 



x 2 > 2 



g -+ f -x 2 



For simplicity we divide into two regimes 



< 



/, 2 
— | X 

9 



> 2 



and 



A, 


> 


/ 




A 9 


1 2 


/ 














< - X 










9 
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9 







> 2. 



Remember that the positivity of flaton masses requires 







/ 


< x 2 








9 



Combining altogether we get 



and 



> 



(63) 

(64) 
(65) 

(66) 

\ in the regime 



> in the regime (| 
fl65l) . Besides, if flatino production rates are sizable, we also have to impose Rp 2 = T(F2 — > 
Xih°)/V(F 2 aF x ) > 10 and M Pl > M Xl + M h o to open the decay mode F x -> Xih°- The 
condition Rp 2 > 10 gives the restriction \f/g\ < 0.02xN Xl and |-| > j^. 



Then we study, in our limit, the constraints given by the conditions R 



h°A )/T{^' ^aF 1 )> 10 and R F . 



R 



R 



Fi 



R 



F 2 




^ = T(tp' - 
> 10. The ratios R are 



(67) 
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|^|xV and M| 2 ~ 12 ^x 2 fi 2 for ^ < 12j|j, and M Fl <-> A7 Fa 



where M| 

12|*|. We can now divide the parameter space into four regions. 



for ^ > 



I) 
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L 


A 4 




9 



X 

■Vi 



, II) 2 < 



< 



Depending on ^ < 12 



,7/7) 



< 



< 2 



< 



we will have region a or region b. 



or — > 12 

M 1 _ 

We find that all of the a regions are forbidden, and so is the 714 region. The constraints 



for the regions are as follows. 

I b x> 14, 

1 < 



Ag < 0, 



12 
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x 2 
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II b x>9, A g > 0, 



2 < 



< 3 x 1(T 2 



Af 




f 




2 < 






f_ 


< 




x\ 




< 








9 






fi 




9 



X 



III b A g >0, 1< 2 



< 



< 2 x 10" 
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< 








9 




fi 



< 2. 



To summarize, we find that: 

- In all cases, x has to be large (> 10) 

- In cases 77 and 777, ^ 
case 7. In all cases x 2 ^ 

9 

- In all cases, \Af \ ~ |A £ 



(68) 



has to be small (< 3 x 10 2 ) but it has no useful upper bound in 
> 1. 

~ \fi | is a possibility. 



7 Conclusions 

We have explored the cosmology of a particularly attractive extension of the Standard Model, 
which has a Peccei-Quinn symmetry broken only by two 'flaton' fields <f)p and 0q, character- 
ized by a very large vev (10 10-12 GeV) and a relatively small mass (10 2 3 GeV). These and 
their superpartners generalize the saxion and axino, that in the non-flat case are the only 
fields with soft masses. 

In contrast with more general models the density of dark matter axions can be estimated 
with essentially no assumption about other sectors of the theory, assuming only that 0p has 
a positive effective mass-squared in the early Universe. If 0q also has a positive effective 
mass-squared the axion is an excellent dark matter candidate. In the opposite case the axion 
density is probably to low in this particular model. 

Our main concern has been with a different, highly relativistic, axion population that is 
produced by flaton decay. We have calculated the rates for all relevant channels and examined 
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the constraint that the energy density of these axions does not upset the predictions of the 
standard nucleosynthesis. 

For the KSVZ model we find that it is almost impossible to satisfy such a cosmological 
bounds due to the fact that the only flatonic decay channel in competition with the axionic 
one is into gluons (through axial anomaly) which is too much constrained. 

For the DFSZ model there are more decay channels. To evade complicated phase space 
suppressions we concentrate on the Higgs decay products, as the mass of the lightest Higgs 
has naturally a relatively low upper bound and the mass of the other Higgs and flaton fields 
are fixed by the parameters of the flatonic potential itself. We have quantified a portion of 
the parameter space available showing the strength of this model. 

An interesting question, lying beyond the present investigation, is whether the allowed 
region of parameter space can be achieved in a supergravity model with universal soft pa- 
rameters. 
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